Updated CMB constraints on Dark Matter annihilation cross-sections 
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The injection of secondary particles produced by Dark Matter (DM) annihilation at redshift 
100 ^ z ^ 1000 affects the process of recombination, leaving an imprint on Cosmic Microwave 
Background (CMB) anisotropies. Here we provide a new assessment of the constraints set by CMB 
data on the mass and self-annihilation cross-section of DM particles. Our new analysis includes the 
most recent WMAP (7- year) and ACT data, as well as an improved treatment of the time-dependent 
coupling between the DM annihilation energy with the thermal gas. We show in particular that the 
improved measurement of the polarization signal places already stringent constraints on light DM 
particles, ruling out 'thermal' WIMPs with mass < 10 GeV. 

PACS numbers: 26.35.-l-c, 98.80.Cq, 98.80.Ft 
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Introduction. Precision measurements of the Cos- 
mic Microwave Background (CMB) anisotropy in tem- 
perature and polarization represent powerful tools to 
constrain new physics processes (see e.g. In par- 

ticular, the remarkable agreement between the theoret- 
ical description of the recombination process, occurring 
at Zr ~ 1000, and CMB data, severely constrains new 
sources of ionizing photons, and more in general any devi- 
ation from standard recombination [2[ , as recently shown 
by several groups of authors (see e.g. 

In our previous paper @ (hereafter GIBM09), we stud- 
ied the constraints that this analysis can set on the prop- 
erties of dark matter (DM) particles [7l-[To|. under the 
assumption that standard recombination is modified by 
dark matter annihilation only. 

Here, we present an update of GIBM09, and obtain 
new constraints on the DM particle annihilation cross 
section and mass, based on more recent data (WMAP 7- 
year [ll| and ACT 2008 [13] data), and on a new and 
more accurate parametrization of the coupling of the 
DM-induced energy shower to the thermal gas. 

Our results turn out to be competitive with constraints 
of diverse astrophysical nature, such as radio observation 
of the galaxy, antiprotons, gamma rays from the Galac- 
tic center and Galactic halo [T3l-[l6|. but with respect to 
them, they have the advantage of not being affected by 
large astrophysical uncertainties. In fact, our CMB con- 
straints arise from redshifts in the range 100 ^ z ^ 1000, 
i.e. well before the formation of any sizable gravitation- 
ally bound structure, and they therefore do not depend 
on highly uncertain parameters related to structure for- 
mation, such as halo shape, concentration or minimal 
mass. 

Annihilating DM and Recombination. We 

briefly recall here the effect of energy injection from DM 
annihilation on the recombination history (see GIBM09 
for further details). High-energy particles injected in the 
high-redshift thermal gas by DM annihilation (or decay) 



are typically cooled down to the keV scale by high energy 
processes (see details below) ; once the shower has reached 
this energy scale, the produced secondary particles can i): 
ionize the thermal gas, ii): induce Ly-a excitation of the 
hydrogen and Hi): heat the plasma; the first two modify 
the evolution of the free electron fraction Xe, the third 
affects the temperature of baryons. The rate of energy 
release ^ per unit volume by a relic self-annihilating 
DM particle is given by 
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where (av) is the effective self-annihilation rate and m-^ 
the mass of the DM particle, Udm the DM density pa- 
rameter and pc the critical density of the Universe today; 
the parameter f[z) indicates the fraction of energy which 
is absorbed overall by the gas, under the approximation 
that the energy absorption takes place locally. We note 
that the presence of the brackets in (av) denote a thermal 
average, as appropriate for relativistic particles at decou- 
pling. At the redshifts of interest here {z ~ 1000) the 
relative velocities of DM particles are v ^ vs = lO^^c, 
i.e. in the extreme non-relativistic limit. Though hold- 
ing for s-wave annihilations, (crw)_8 ~ {(Jv)dec, the same 
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is not true in general. For instance, (crw)- 
in models with so-called Sommerfeld enhancement, and 
{av) -R'^,{av)dec in models with p-wave annihilations 

(e.g. Wm- 

In GIBM09, we considered the fraction of energy /(z) 
absorbed by the plasma to be constant with redshift, 
f{z) = f. In the following sections, we will present up- 
dated constraints obtained by supposing that / is con- 
stant with redshift, as well as constraints considering the 
actual DM model dependent redshift shape of f{z), as 
calculated in Slatyer et al (Tgj . 

The formalism we use to introduce the extra energy 
terms in the recombination equations are the same as in 
GIBM09, but here we additionally consider the modifica- 
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tions to both hehum and hydrogen recombination and we 
change the extra ionization term (Eq. 6 in GIBM09). In 
this equation we do not include Peebles' C factor, since 
we assume that the extra ionization photons only ionize 
the ground state. We checked that excluding the factor C 
in this term in any case does not change the constraints 
by more than 6% in the case of a constant /. The ex- 
tra Lyman alpha term remains the same as in Eq. 7 of 
GIBM09. 

The most remarkable effect of injecting energy around 
the recombination epoch is that the amount of free elec- 
trons that survive at low redshift after recombination is 
larger compared to the standard case where no annihila- 
tion happens. The CMB spectra are therefore affected as 
follows (see Ref. [l^ for a more detailed discussion) : the 
enhanced amount of Xe at low redshift increases the width 
of last scattering surface, and consequently the width of 
the visibility function. This results in a suppression of 
the amplitude of the oscillation peaks in the tempera- 
ture and polarization power spectra, especially at scales 
smaller than the width of the last scattering surface. 

This effect is degenerate with that of other cosmo- 
logical parameters, affecting the amplitude of peaks at 
low/high multipoles, such as the scalar spectral index 
Us and the baryon density cjf, = ftbh^. In particular, a 
value of Pann different from zero can be compensated by 
a higher value of the scalar spectral index Ug , that gives 
more power to smaller scales of the spectrum and, to 
a much smaller extent, by a higher value of the baryon 
density Wf,. Nevertheless, this last parameter changes the 
relative amplitudes of the peaks and is therefore less de- 
generate with Pann than Us ■ A smaller degeneracy is also 
found with the dark matter density lUc = ^ch^- 

On the other hand, a larger width of the recombination 
epoch due to DM annihilation increases the quadrupole 
moment of the radiation field as well, enhancing the 
amplitude of the polarization power spectrum at large 
scales. Furthermore, the fractional contribution to the 
quadrupole of the temperature monopole with respect 
to the temperature dipole is increased, therefore slightly 
shifting the position of the peaks of the polarization 
power spectrum (see [l^l for further details). An ac- 
curate measurement of the polarization power spectrum 
can therefore help breaking the degeneracy with n^. 

Updated CMB constraints Following GIBM09, 
we compute here the theoretical angular power in pres- 
ence of DM annihilations, by modifying the RECFAST 
routine [2l| , along the lines described in the previous sec- 
tion, in the C AMB code [l^] , and by making use of pack- 
age cosmomc ^23 ] for Monte-Carlo parameter estimation 
(see GIBM09 for the details of our statistical analysis, 
including convergence tests and priors). The fractions of 
energy release inducing ionization, Lyman-a excitation 
or heating of the baryonic gas are the same as in Chen 
and Kamionkowski's pape r IT7| . based on the results of 
Shull and Van Steenberg [l^. The dependence on the 
properties of the DM particles is encoded in the quantity 
Pann, appearing in eq. [U that we use as a parameter 



in the code. We start by considering / as a constant pa- 
rameter, in order to quantify the impact of the new CMB 
data sets on the analysis presented in GIBM09. 

Besides Pann 1 we sample the following six-dimensional 
set of cosmological parameters (with flat priors): the 
physical baryon and CDM densities, Wfc and cjc, the scalar 
spectral index, n^, the normalization, hilO^'^As{k = 
0.002/Mpc), the optical depth to reionization, r, and the 
Hubble parameter Hq. We consider purely adiabatic ini- 
tial conditions. 

We include the seven-year WMAP data (WMAP7) 
prf (temperature and polarization), and the 2008 ACT 
telescope data (temperature only), that probes the 
temperature angular power spectrum at small scales 
(600 I ^ 8000, but in our analysis we consider 
^max = 3500). The routines for computing the likeli- 
hoods for each experiment are supplied by the WMAP 
and ACT teams respectively. 

The results of our analysis are in TableU there we show 
the constraints on the DM annihilation parameter Pann 
and on the parameters that are more degenerate with 
it (i.e. ns,uJb,ujc) obtained using WMAP7 data and the 
combination of WMAP7 plus ACT data. We also report 
the the constraints obtained by the WMAP7 team [23| 
and by the ACT team 25] on the cosmological parame- 
ters in the standard case (no DM annihilation) in order to 
show the bias introduced on the cosmological parameters 
by not considering DM annihilation in the analysis. 

Note in particular that the constraint on Pann with 
WMAP7 data is improved by a factor ^1.8 with respect 
to the WMAP5 constraint obtained in GIBM09. The 
7-year data release has in fact a better measurement of 
the third peak of the temperature power spectrum at 
I ^ 1000—1200 and of the second dip in the temperature- 
polarization power spectrum at I ^ 450. This allows a 
better measurement of ujb and Uc and a partial break of 
the degeneracy with Pann- On the other hand, the bias on 
Us remains noticeable at 1 — cr level, as the measurement 
of the polarization power spectrum is still not sufficient to 
break the degeneracy with Pann ■ Adding the information 
at small scales from the ACT data additionally improves 
the constraint on Pann by 13%. As we can see from the 
table, the Harrison-Zcl'dovic model n = 1 is consistent 
with the WMAP+ACT analysis in between two standard 
deviations when dark matter annihilation is considered. 

We have also checked our results when other non- 
standard parameters are considered in addition to the 
standard ones and Pann- We have considered one of the 
following additional parameters at the time: the fraction 
of Helium abundance Yhe , the massive neutrino density, 
ri^, and the running of the scalar spectral index a. None 
of these parameters appeared to be degenerate with Panm 
therefore not affecting the results on the upper limits re- 
ported in Table m 

Implementing the redshift dependence of /. 
We have so far worked under the assumption that the 
fraction of the rest DM mass energy absorbed by the 
plasma is constant with redshift. Yet, the fraction of 
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WMAP7 


WMAP7+ACT 


WMAP7 Standard WMAP7±ACT Standard 


Pann[cm^ /s/GeV] 


< 2.42 X 10"^^ 


< 2.09 X lO"^'^ 




Us 


0.977 ±0.015 


0.971 ±0.014 


0.963 ±0.014 0.962 ±0.013 


lOOfib/i^ 


2.266 ± 0.057 


2.237 ± 0.053 


2.258lo;j^6 2.214 ±0.050 




0.1115 ± 0.0054 


0.1119 ±0.0053 


0.1109 ± 0.0056 0.1127 ±0.0054 



TABLE I: Constraints on the annihilation parameter pann and on the cosmological parameters that are more degenerate with 
it, i. e. the scalar spectral index Us, the baryon density ujb and the dark matter density ujc- We report the results using WMAP7 
data and WMAP7±ACT data. The constraints on pann are upper bound at 95% c.L, while for the other parameters we show 
the marginalized value and their errors at 68% c.l. The last two columns reports the value of the cosmological parameters in 
the standard ACDM case with no annihilation, as found by the WMAP7 team [ij and the ACT team [25l |. 



the initial energy deposited into the gas is not constant 
with cosmic time, even if the on-the-spot approximation 
holds true at all redshifts of interest. This problem has 
been addressed in [T^ , where the authors have computed 
the evolution of the energy fraction f{z) for different pri- 
mary species, and DM particle mass. As it can be seen 
from their Figure 4, the f{z) is a smoothly varying func- 
tion of rcdshift (even more so for the values of interest in 
our problem 100 % z 1000). We show the constraints 
for time- varying f{z) in Figure [1] Interestingly, the new 
results rule out 'thermal' WIMPs with mass < 10 
GeV. 

We have checked the constraints which is possible to 
place using the redshift dependent shape of / presented in 
Equation Al and Table 1 of [l^. We have obtained con- 
straints for purely DM models annihilating solely (and 
separately) into electrons and muons, with different DM 
masses, reported in Table [TTl This choice of annihila- 
tion channels brackets the possible values of f{z): the 
case of annihilation to other channels (except of course 
neutrinos, which practically do not couple at all with the 
plasma) falls between the two limiting cases studied here. 

Although the implementation of the z-dependence of 
/ clearly leads to more accurate results, we found that 
taking a simplified analysis with constant /, such that 
f{z — 600) — f const, leads to a difference with respect to 
the full /(z) approach of less than ^ 15%, depending on 
the annihilation channel considered. 

Discussion and Conclusions. In this brief report 
we have provided new updated CMB constraints on 
WIMP annihilations, with an improved analysis that 
includes more recent CMB data (WMAP7 and the 
ACT2008) and implementing the redshift evolution of the 
thermal gas opacity to the high energy primary shower. 
We have also found that a simplified analysis with con- 
stant / = f(z — 600) leads to an error on the maximum 
DM self-annihilation cross section smaller than ~ 15%, 
with respect to a treatment that fully takes into account 
the redshift dependence of f{z). 

While we were finalizing this paper, Hutsi et al. 
(IICIIR2011) [1^ have reported results from a similar 
analysis, using an averaged evolution of the f{z). They 
provide 2 — a upper limits from WMAP7 with 1 — a un- 
certainties on these limits due to the method used. These 
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FIG. 1: Constraints on the cross section < av > in function 
of the mass, obtained using a variable f{z) for particles anni- 
hilating in muons {x signs) and in electrons (diamonds) using 
WMAP7 data (red) and WMAP7±ACT data (black) at 95% 
c.L. The exclusion shaded areas are obtained for interpolation 
of the WMAP7 ± ACT data points for muons (dark shading) 
and electrons (light shading). The black solid line indicates 
the standard thermal cross-section < av >= 3 x lO"'^^ cm^/s. 



results are a factor between 1.2 and 2 weaker than ours. 

This is partially due to the fact that we account for ex- 
tra Lyman radiation in our code, but this can account for 
only less than 10% of the difference between the results. 

As in GIBM09, we have calculated how much the 
Planck satellite and a hypothetical Cosmic Variance Lim- 
ited experiment will improve the constraints compared 
to WMAP7 in the case of constant / (constraints for 
Planck and CVL reported in GIBM09). We obtain im- 
provement factors of 8 and 23 for Planck and CVL re- 
spectively, which are compatible with the ones reported 
in HCHR2011, 6 and 13. The difference for the CVL 
experiment is attributed to the slightly different specifi- 
cations used for the CVL experiment in IICHR2011 and 
in GIBM09, namely the maximum multipole considered 
in the analysis, as also stated in HCHR2011. Clearly the 
data from the on-going Planck satellite mission, expected 
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< av > in [cm''/s] with Variable / 



< av > in [cm'*/s] with Constant / = f{z = 600) 



mx[GeV] channel 


WMAP7 




WMAP7+ACT 


f{z = 600) 


WMAP7 




WMAP7+ACT 


1 GeV 


e"^e~ 


< 2.90 X 10" 


27 


< 2.41 X 10'^^ 


0.87 


< 2.78 X 10" 


27 


< 2.41 X lO"^'^ 


100 GeV 


e'''e~ 


< 3.95 X 10" 


25 


< 3.55 X 10"^^ 


0.63 


< 3.87 X 10" 


25 


< 3.35 X lO"^'^ 


ITeV 


e'''e~ 


< 4.68 X 10" 


24 


< 3.80 X 10"^* 


0.60 


< 4.02 X 10" 


24 


< 3.48 X 10"^" 


1 GeV 




< 8.68 X 10" 


27 


< 6.93 X 10"^^ 


0.30 


< 8.03 X 10" 


27 


< 6.95 X lO-^'^ 


100 GeV 




< 9.82 X 10" 


25 


< 8.94 X 10"^^ 


0.23 


< 1.03 X 10" 


24 


< 8.91 X 10"^^ 


ITeV 


11+ ii" 


< 1.20 X 10" 


23 


< 9.41 X 10"^"* 


0.21 


< 1.15 X 10" 


23 


< 9.96 X 10"-^'' 



TABLE II: Upper limits on self-annihilation cross section at 95% c.l. using WMAP7 data and a combination of WMAP7 and 
ACT data . On the left-side of the table, we show the results obtained using the proper variable f{z) for each model. On the 
right side, for sake of comparison, we show the results obtained by taking the constraints for a constant generic / reported in 
Table [H and then calculating < gv > for each case imposing that / is equal to the corresponding f{z — 600) for each model. 
We show results for particles annihilating in electrons and muons. 



to be released by early 2013, will play a crucial role in 
constraining additional sources of ionization, such as DM 
annihilation, in the early universe. 
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